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ABSTRACT 

Measurements of the two-point angular correlation function w(6) for 407 short gamma- 
ray bursts collected in the Current BATSE Catalogue reveal a ~2<r deviation from 
isotropy on angular scales 8 ~ 2 — 4 degrees. Such an anisotropy is not observed in the 
distribution of long gamma-ray bursts and hints to the presence of repeated bursts for 
up to ~ 13% of the sources under exam. However, the available data cannot exclude 
the signal as due to the presence of large-scale structure. Under this assumption, 
the amplitude of the observed w{ff) is compatible with those derived for different 
populations of galaxies up to redshifts ~ 0.5, result that suggests short gamma-ray 
bursts to be relatively local sources. 
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The Burst And Transient Source Experiment (BATSE) on 
board the Compton Gamma Ray Observatory collected 
more than 2000 Gamma-Ray Bursts (GRBs) in its 9-year 
activity. The analysis of such a large sample has deeply 
characterized the temporal and spectral properties of GRBs 
(Fishman & Meegan 1995). Despite the multitude and com- 
plexity of the light curves observed, the burst duration, typ- 
ically defined by T90 as the time interval during which 90% 
of the fluence is accumulated (Kouvelitou et al. 1993), in- 
dicates the existence of two main classes: short (T90 < 2 
s) and long (T90 < 2 s) bursts. This distinction is also rein- 
forced by the associated spectral shape since short bursts are 
typically harder than long ones (Fishman & Meegan 1995; 
Mitrofanov et al. 1998). At present, it is not known if this 
dichotomy is produced by different progenitors or by differ- 
ent physical mechanisms originating from the same object 
or, again, from different observing conditions. Nonetheless, 
it has been suggested (MacFadyen & Woosley 1999) that 
long bursts might be produced by the core collapse of a 
single massive stars (the Collapsar model - Woosley 1993; 
Paczynski 1998) whereas short bursts might be more effi- 
ciently generated by double compact object mergers (Good- 
man 1997; Narayan, Piran & Kumar 2001) or double neutron 
star heating (Salmonson & Wilson 2002). A distinctive fea- 
ture of these two scenarios is the expected GRB distance: 
compact object mergers should have a lower median redshift 
than Collapsars (Fryer, Woosley & Hartmann 1999) which 
would be tracers of the high redshift universe (Lamb 2001). 



A cosmological origin of long bursts was first suggested 
by the isotropic sky distribution observed by BATSE (Briggs 
1993; Meegan et al. 1992; Paczynski 1995) but the definitive 
proof came with the first redshift measurement for GRB 
970508 (Metzger et al. 1997) at z = 0.835. At present, the 
distribution of known redshifts extends up to z=4.5 with a 
mean (z) ~ 1.27 (Andersen et al. 2000; Djorgovski et al. 
2001). However very little is presently known about short 
GRBs because no radio, optical or X-ray counterpart has 
been associated with any burst lasting less than 2 s (Hurley 
et al. 2002), although the possible existence of short burst af- 
terglows has been suggested (Lazzati, Ramirez- Ruiz & Ghis- 
ellini 2001). Thus, the distance scale of short GRBs is still 
an open issue and its solution would have important impli- 
cations for the understanding of their progenitors. 

The aim of this Letter is to investigate the nature of 
short GRBs by means of statistical techniques such as the 
angular two-point correlation function, since different hy- 
potheses are expected to leave different imprints on their 
distribution on the sky. In particular we will concentrate on 
the two distinct possibilities of burst repetition and large- 
scale structure. 

In order to do so, we have considered sources from the public 
BATSE on line catalogue* which contains 2702 bursts from 
GRB 910421 to GRB 000526. For our analysis, a key ingre- 
dient is the estimate of the duration T90, which has been 
derived only for 2041 events of this catalogue. Also, the flux 
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and fluence are only available for 2135 events. The maxi- 
mum number of bursts that we could extract with measured 
durations and fluences is 2035, of which 497 are short events 
with T 9 o < 2 s. 

Our work is presented as follows: in Section 2 we briefly 
recall the results of previous works examining evidence for 
GRB anisotropics, while in Section 3 we introduce the 
method and show the results on the angular correlation 
function. Section 4 discusses some possible models for the 
observed distribution of short GRBs and summarizes our 
conclusions. 



2 EVIDENCE OF GRB ANISOTROPIES 

First studies on the angular distribution of GRBs by Hart- 
mann & Blumenthal (1989) including 80 GRBs detected by 
the IPN and 160 GRBs collected by KONUS, showed a pos- 
sible deviation from isotropy on scales ~ 5°. A more re- 
fined analysis (Quashnock & Lamb 1993) based on the first 
BATSE (IB) catalogue (260 bursts) confirmed this evidence 
but with a cut on the burst positional errors A < 10° (but 
see cautions in Narayan & Piran 1993; Brainerd 1996). 

The analysis of joint spatial and temporal clustering of 
the bursts in the IB catalog (Wang & Lingenfelter 1995) 
showed an excess of bursts with angular separation lower 
than 4° within 4 days with a small probability (1 x 10~ J ) 
of random occurrence. However, less significant results are 
claimed by Petrosian & Efron (1995). In particular, the need 
of larger samples for investigation of the possible burst rep- 
etition by means of the two-point correlation function is 
discussed by Brainerd (1996). 

Evidence of possible burst clustering was recently 
claimed for both the populations of short and intermedi- 
ate GRBs from the analysis of larger samples (Balazs et al. 
1998; Meszaros et al. 2000). The possible different nature 
of the short and long events, in fact suggests the separate 
analysis and comparison of their clustering properties. 



3 THE METHOD 

One of the possible statistical tests for the detection of 
anisotropy is provided by the two-point angular correla- 
tion function w(Q), which gives the excess probability, with 
respect to a random Poisson distribution, of finding two 
sources in the solid angles 5Q.i 5Q.2 separated by an angle 
9. We applied this technique to the distinct samples of long 
and short GRBs in order to investigate for possible devia- 
tions from isotropy in their spatial distribution. 

Given the large uncertainties affecting the determina- 
tion of the positions of GRBs, especially in the case of short- 
lived ones (see Figure 1), we have decided to only consider 
those objects which have angular coordinates with associ- 
ated errors A < 10°. The choice for this value results as 
the best compromise between the requirements of precision- 
position measurements and large-enough datasets to allow 
for sensible statistical analysis. By applying the above cut, 
we end up respectively with N\ oa& = 1486 (out of 1538) and 
-Wahort — 407 (out of 497) sources in the case of long-lived 
and short-lived GRBs. 

Random catalogues have subsequently been generated 
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Figure 1. Distribution of positional errors for objects included 
in the Current BATSE Catalogue. The dotted line is for the pop- 
ulation of long-lived GRBs, while the solid one is for short-lived 
bursts. The vertical line indicates the upper limit for inclusion in 
the datasets considered in this work (see text for details). 



for each of the two subsets, with positions of the objects 
uniformly distributed on the celestial sphere. In order to 
minimize Poisson fluctuations due to the limited number 
of sources in the original datasets, random catalogues were 
required to include 10 x Ai on g objects in the case of long- 
lived GRBs and 20 x A^ort objects in the case of short-lived 
GRBs. We then counted the number of distinct data-data 
pairs (DD), data-random pairs (DR), and distinct random- 
random pairs (RR) as a function of angular separation, and 
calculated the two-point angular correlation function w(9) 
using the estimator (Hamilton 1993) 

«(*) = a -jW - 1 « 

on the angular scales 1 < 9 < 20 degrees. 

The open dots in Figs. 2 and 3 illustrate the results re- 
spectively for the population of long-lived and short-lived 
GRBs; the error bars indicate Poisson estimates for the 
points. We note that, while in the first case values of w(9) 
at all angular scales seem to be consistent with an isotropic 
distribution, the angular correlation function of short-lived 
GRBs shows a significant excess of power on scales ~ 2 — 4 
degrees. 

In order to test the significance of this result, we have 
furtherly examined two reduced datasets, obtained from the 
original one by only considering those sources with posi- 
tions known to better than A < 5° (1239 long and 196 
short GRBs). Solid squares in both Figs. 2 and 3 illustrate 
the values of w(9) obtained in this case. Once again, mea- 
surements of the angular correlation function for long-lived 
GRBs are entirely consistent with absence of intrinsic non- 
randomness, in agreement with earlier results (e.g. Meegan 
et al. 1995, 1996; Bennett & Rhie 1996; Tegmark et al. 1996) 
obtained for samples which mainly included this population 
of gamma-ray emitters, and with more recent findings (e.g. 
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Figure 2. Angular correlation function w(0) of long-lived 
gamma-ray bursts from the Current BATSE Catalogue. Empty 
dots are for sources with positional uncertainties A < 10° (1486 
objects), while filled squares are for sources with positional un- 
certainties A < 5° (1239 objects). 
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Figure 3. Angular correlation function w(6) of short-lived 
gamma-ray bursts from the Current BATSE Catalogue. Empty 
dots are for sources with positional uncertainties A < 10° (407 
objects), while filled squares are for sources with positional un- 
certainties A < 5° (196 objects). 



Balazs et al. 1998; Meszaros et al. 2000) concerning the pop- 
ulation of long-lived GRBs alone. 

On the other hand, the distribution of short-lived, A < 
5° GRBs - even though with larger errors due to the small 
number of objects considered in the analysis - shows a de- 
parture from isotropy on small angular scales even stronger 
than it was for the previously-analysed A < 10° subsample. 
We can then conclude that the distribution of short-lived 
gamma-ray bursts exhibits a ~2u deviation from isotropy 
on scales ~ 2 — 4 degrees. 



4 DISCUSSION 

Given the above result on the anisotropic distribution of 
short-lived GRBs, it is worthwhile investigating (and pos- 
sibly excluding) all the systematics which could produce a 
non-zero signal in the angular correlation function. 
One of the possible causes for systematic deviations from 
isotropy is related to the non-uniform sky-exposure of 
BATSE (Fishman et al. 1994; Band 1996), which in prin- 
ciple could lead to measurements of anisotropy even for ran- 
domly distributed sources. We note however that such an 
effect would affect the observed distribution of both long 
and short GRBs in a similar fashion, since there is no a pri- 
ori reason why different sub-classes should be observed by 
the instrument in different ways. The fact that long-lived 
sources exhibit a w(9) which is highly consistent with the 
hypothesis of isotropy at all the angular scales probed by 
our analysis then tends to exclude responsibility from in- 
strumental effects on the observed anisotropy of short-lived 
GRBs (see also Balazs et al. 1998; Meszaros et al. 2000). 

Another possible cause for the observed anisotropic dis- 
tribution can be attributed to intrinsic Poissonian fluctua- 
tions of short-lived GRBs due to their small number. Such 
an effect would however produce random fluctuations of the 
measured correlation function at all angular scales, while the 
data in our case show that an appreciable deviation from 
zero values is only obtained on small distances; measure- 
ments of w(9) at all 9 > 5° are remarkably consistent with 
a null clustering signal. As a further test, we have re-binned 
sources belonging to the sub-class of short-lived GRBs and 
evaluated the two-point correlation function by means of 
equation (1) in the angular range 1 < 9 < 10 degrees. 
Results from this analysis are presented in Figure 4: once 
again a marked deviation from isotropy is only detected on 
scales 9 ~ 2 — 4 degrees, in full agreement with the findings 
of the previous Section. 

Having excluded the effects of systematics on the de- 
tected anisotropy, we are then left with intrinsic non- 
randomness in the distribution of short-lived GRBs. This 
intrinsic anisotropic distribution would arise if e.g. some of 
the sources were repeating bursts, as indeed expected in sev- 
eral GRB models (Quashnock & Lamb 1993; Wang & Lin- 
genfelter 1995; Meegan et al. 1995; Bennett & Rie 1996; 
Tegmark et al. 1996). 

If we define / as the fraction of observed bursts that 
can be labeled as repeaters and v as the average number 
of observed events per repeating source (y > 2) then, for 
positional errors drawn from a Gaussian distribution with 
standard deviation £ = A/(2V21n2) << 1 (A ~ 5°, see 
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Figure 4. Angular correlation function w(0) for short-lived 
gamma-ray bursts with positional uncertainties A < 10°. The 
dashed line shows the best fit to the data obtained under the hy- 
pothesis of burst-repetition, while the dotted line illustrates the 
best case for large-scale clustering (see text for details). 



Figure 1), the correlation signal arising from the repeating- 
burst assumption can be written as (see Meegan et al. 1995) 



w(6) 
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Equation (2) was used to fit the data presented in Figure 
4 for different values of the product / (y — 1). It turns out 
that the best-fit value is given by f(v— 1) = O.OTIq^, result 
which rejects the hypothesis of non-repetitivity at the ~ 2a 
confidence level. 

By then fixing the average number of repeating events, we 
obtain the following constraints on the fraction of sources 
expected to repeat: 



o.o7lg:£ 



v = 2^ f 
i/ = 3 - / = 0.035tH? 
v = 4 - / = 0.023^3 

In other words, the data indicate that up to 13 per cent 
of short GRBs may have burst on average twice during the 
years corresponding to the lifetime of BATSE. This fraction 
obviously decreases for an increasing number of repetitions, 
going down to a more negligible ~ 2 per cent in the case of 
v = 4. 

The above explanation is however not the only one ca- 
pable to account for the observed anisotropy on a few de- 
grees scale. Another alternative which needs to be taken into 
account is the presence of large-scale structure, arising un- 
der the hypothesis that short GRBs trace the distribution of 
the galaxies in which they reside (note that here we are ex- 
plicitly assuming GRBs to be extragalactic sources). In this 
framework, the angular correlation function can be modeled 
(to a first approximation) as the power-law w(9) = A O 1 - 1 ^^ 
(see e.g. Peebles 1980). 

Fits to the data presented in Figure 4 give values for the 



amplitude of the correlation function A ~ 0.3±0.2, indepen- 
dent (within the errors) of the different values for the slope 
7 = 1.6,2.0,1.8, chosen so to mimic the clustering prop- 
erties respectively of star-forming and early-type galaxies, 
and of a mixture of both. The case for 7 = 1.8 is illustrated 
in Figure 4 by the dotted line. Despite the large errors as- 
sociated to its determination, the allowed range of values 
for A is in good agreement with results e.g. from the APM 
Survey (Loveday et al. 1995) based on samples of 2 < 0.2 
galaxies divided by morphological type. On the other hand, 
clustering analyses performed at higher redshifts find am- 
plitudes which are much smaller than the values indicated 
by the GRB data (A < 10~ 2 for z ~ 1 - 1.5 EROs, Roche 



et al. 2002; A 



10~ J for 2 



3 Lyman Break galaxies, 



Porciani & Giavalisco 2002; A ~ 5 • 10^ for 2 ~ 2.5 - 3.5 
galaxies found in the Hubble Deep Field, Magliocchetti & 
Maddox 1999). If we then use large-scale structure to ex- 
plain the clustering measurements presented in this work, 
we can conclude that - if confirmed - they imply short-lived 
GRBs to be hosted by relatively (possibly 2 < 0.5) local 
galaxies, in agreement with the median redshift 0.5-0.8 pre- 
dicted by neutron star/neutron star or black hole/neutron 
star merger models (Fryer et al. 1999). In fact (see e.g. Mac- 
Fadyen & Woosley 1999), the alternative to produce short 
bursts in the context of the Collapsar model would meet 
some difficulties at 2 < 0.8 since the intense mass loss from 
the progenitor - due to a higher metallicity with respect to 
that of stars at higher redshifts - could prevent the collapse 
of the helium core. 

Obviously, one would like to discern between an 
anisotropic signal caused by burst repetition and one due to 
large-scale structure. Unfortunately, the still small number 
of detected sources does not allow to probe angular scales 
9 < 1 degrees where these two effects would greatly dif- 
ferentiate from one another. The solution to this problem is 
however within reach, since the goal of collecting hundreds of 
GRBs per year will be soon achieved by SWIFT (Burrows et 
al. 2002) , which will also allow for precise position measure- 
ments in order to perform optical to X-ray spectrophotome- 
try, with the aim to determine the origin of GRBs and their 
afterglows and to use bursts to investigate the early Uni- 
verse. Moreover, the likely detection of several short bursts 
(which was severely limited in the BATSE experiment due to 
its lower sensitivity to extremely fast transients) by INTE- 
GRAL and their location with sufficient accuracy through 
the IPN system will lead to a better understanding of the 
still mysterious population of short bursts. 
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